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OVERALL GOALS 


The introduction of composite materials is having a profound effect on aircraft design. Since 
these materials permit the designer to tailor material properties to improve structural, aerody- 
namic and acoustic performance, they require an integrated multidisciplinary design process. Fur- 
thermore, because of the complexity of the design process numerical optimization methods are 
required. 

The utilization of integrated multidisciplinary design procedures for improving aircraft design 
is not currently feasible because of software coordination problems and the enormous computa- 
tional burden. Even with the expected rapid growth of supercomputers and parallel architectures, 
these tasks will not be practical without the development of efficient methods for cross-disciplinary 
sensitivities and efficient optimization procedures. 

The present research is part of an on-going effort which is focused on the processes of simul- 
taneous aerodynamic and structural wing design as a prototype for design integration. A sequence 
of integrated wing design procedures has been developed in order to investigate various aspects of 
the design process. 


o NEED 

• composite materials 

• aircraft complexity 

o PAY-OFF 

• better designs 

O DIFFICULTIES 

• computational cost 

• software coordination 
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PREVIOUS RESEARCH EFFORTS 


In their initial efforts, the authors considered the integrated design of a high aspect-ratio 
sailplane wing. The sailplane mission was used to illustrate the advantages of including aerodynamic 
and structural interactions in the design process, by optimizing for circling flight in a thermal 
current followed by cross-country cruise. Furthermore, the simplicity of the sailplane wing planform 
and structural design allowed for the use of rudimentary analysis methods, (lifting-line and beam 
theory). The simplicity of these analyses made it feasible to calculate all the sensitivity derivatives of 
the aerodynamic shape and structural sizes, along with all the cross-sensitivity derivatives, directly, 
without any further approximation, at each step of the numerically optimized design process. The 
results, reported in Ref. 1, demonstrated that integrating the structural and aerodynamic design 
processes leads to wing designs superior to those obtained by the traditional sequential approach. 

The next step of the integrated wing design procedure study again involved the sailplane wing 
design, but with analysis methods which are representative of methods used for low-speed aircraft 
wing designs. The utilization of a vortex-lattice method and a structural finite-element method, 
while providing for a more exact analysis and allowing for more general wing shapes, introduced 
the need for more design variables and constraints, and were significantly more expensive to use in 
the design process. In Ref. 2, it was shown that by incorporating perturbation methods for cross- 
sensitivity calculations and approximate optimization procedures, an estimated 10 hours of IBM 
3084 CPU time for a complete integrated design, was reduced to less than ten minutes. Most of the 
remaining computational cost was associated with the calculation of derivatives of the aerodynamic 
influence coefficient matrix and the structural flexibility matrix. 


O Demonstrated benefits of integrated design using rudimentary analysis meth- 
ods for sailplane design 

O Reduced computational costs by approximate optimization 

O Computational costs remain high due to sensitivity derivatives of aerodynamic 
and flexibility matrices 
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PRESENT OBJECTIVE 


The present paper represents the third step of this study. The objective here is to develop 
an integrated wing design procedure for a subsonic transport aircraft. We still use vortex-lattice 
aerodynamics (so that we are restricted to subsonic speeds) and finite-element structural analy- 
sis. Even with basic aerodynamic design variables, (planform shape and twist distribution), the 
increased complexity of an integrated transport design over the previous sailplane wing design 
requires further computational reductions. We consider two approaches for reducing the computa- 
tional burden of multidisciplinary optimization: 

i. the development of efficient methods for cross-sensitivity calculation; and 

ii. the use of approximate optimization procedures. 

The sensitivity calculation is based on a modular sensitivity method (Ref. 3) for computing 
sensitivity derivatives of a system via partial derivatives of the output with respect to input and 
to design variables of each component of the system. This modular approach, corresponds to the 
abstraction of a system as an assembly of interacting black boxes. This method was developed 
for calculating system sensitivity without modifying disciplinary black-box software packages, Ref. 
4. It allows for the calculation of sensitivity derivatives of a system with a higher accuracy and, 
in most cases, at a lower cost than with conventional finite differencing. The system sensitivity 
derivatives may be used to guide a formal optimization and a Newton’s method solution of the 
coupled interdisciplinary equations describing the system behavior. Within this framework, we 
show that the sensitivities can be computed without the expensive calculation of the derivatives of 
the aerodynamic influence coefficient matrix, and the derivatives of the structural flexibility matrix. 

Furthermore, the same process enables the determination of the sensitivity of the aeroelastic 
divergence dynamic pressure without the determination of the derivatives of the aerodynamic influ- 
ence coefficient matrix and flexibility matrix. This feature should be useful, not only in problems 
of complete integrated aircraft design, but also in aeroelastic tailoring applications. 


O Develop an integrated design procedure for a transport aircraft 
O Utilize a modular sensitivity analysis 
O Reduce computational costs 
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WING DESIGN VARIABLES 


We consider the optimum design of an aircraft wing. The objective function can be the 
structural weight of the wing, an aerodynamic performance index such as the lift-to-drag ratio, 
L/D or a combination thereof. In the present study we minimize the structural weight of the 
wing. The design variables associated with the aerodynamic design include the planform shape 
parameters defined on the figure below, and the twist schedule along the span. 

For the present, preliminary study of integrated structural-aerodynamic design, we assume 
the airfoil shape to be supplied, along with known section characteristics. The design variables 
associated with the structural design are the structural sizes including panel thicknesses and spar- 
cap cross-sectional areas. The finite-element model of the wing is shown schematically in the figure 
below. Additionally, composite material ply orientations in the cover panels are used as design 
variables. 
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INTEGRATED DESIGN PROBLEM 


Constraints are placed on the magnitudes of stresses and strains in the structure, on the aeroe- 
lastic divergence speed, and on aerodynamic performance measures and stall conditions. Additional 
geometric constraints are imposed on the planform shape design variables to prevent unreasonable 
geometries. 

The aerodynamic and structural response is calculated from a coupled set of equations dis- 
cussed below. Aerodynamic performance is calculated at the cruise condition, while the limits on 
stresses and strain are applied for a high-g pull-up maneuver. 


o OBJECTIVE FUNCTION 

minimize wing weight 

o DESIGN VARIABLES 

• planform shape parameters 

• panel thickness , spar cap areas 

• ply orientation 

O CONSTRAINTS 

• stress, strain 

• divergence speed 

• performance, planform shape 
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AEROELASTIC FORMULATION 


The aeroelastic analysis of the wing is simplified by making several assumptions. We assume 
that the effect of the aerodynamics on structural deformations can be approximated by lumping 
the aerodynamic forces at n/ structural grid points (called here the load set), and including only 
the vertical components of the loads. The vector of vertical aerodynamic loads is denoted as F a . 
We assume that the overall aircraft response affects the wing only through the root angle of attack 
a. Finally, we assume that the effect of structural deformations on the aerodynamic response can 
be approximated in terms of the vector of vertical displacements 0 at the load set. 

The vertical aerodynamic loads at the load set, F a , are determined from an aerodynamic 
analysis procedure. For low speed wing designs, we utilize a vortex lattice method (e.g., Ref. 5) 
to compute the lift and the induced drag. The wing is discretized into panels, with each panel 
containing an element of a horseshoe vortex of strength 7 y. By enforcing flow tangency at each 
panel, a vector of circulation strengths T may be computed from eq. (1), below, where p is a 
vector of design parameters and V is a matrix of influence coefficients. The aerodynamic forces are 
computed from a local application of the Kutta-Joukowski theorem, and compressibility effects are 
included through a Gothert transformation. The profile drag for each wing section is calculated 
from the measured airfoil drag polar. The load vector F a is then obtained as eq. (2). Altogether we 
combine equations (l) and (2) as eq. (3), below. The angle of attack is obtained from the overall 
vertical equilibrium of the aircraft as eq. (4), where N is a summation vector, n is the load factor 
and W is the weight of the aircraft. 


• aerodynamic equations 


V(p,6)Y = C(p,ot,6) ( 1 ) 

F a = F a (p,ot,0,Y) ( 2 ) 

Fa = fi{p,<*,0) ( 3 ) 

• vertical equilibrium equation 

h (P, K) = \nW - N t F„ = 0 ( 4 ) 
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STRUCTURAL ANALYSIS 


The vertical displacements at the load set are calculated by finite-element analysis using a 
modification of the WIDOWAC program (Ref. 6). First the nodal displacement vector U is 
calculated by solving eq. (5), where K is the stiffness matrix, T is a Boolean matrix which expands 
F a to the full set of structural degrees of freedom, and Fj is the gravitational and inertia load 
vector. Strains and stresses are then calculated from the displacement vector U. The vertical 
displacements at the load set 0 are extracted from U as eq. (6). Equations (5) and (6) can be 
combined as eq. (7). 


• structural equations 


K(p)U = TF a + nF r (p) ( 5 ) 

0 = T T U ( 6 ) 

8 = h{p,F a ) ( 7 ) 
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SOLUTION PROCEDURE 


Equations (3), (4) and (7) are a set of nonlinear coupled equations for the vector of vertical 
aerodynamic loads, F a , the wing root angle of attack, a and the vector of vertical displacements, 
0. For the analysis problem, the vector of design parameters, p, is given. Reference 3 presented 
a modular sensitivity analysis of such coupled interdisciplinary equations. The modular approach 
permits treating the individual discipline analysis procedures as black boxes that do not need to 
be changed in the integration procedure. Here we employ a similar approach for the sensitivity 
analysis below, with f\ representing an aerodynamic black box and fz a structural black box. We 
also use the same approach for the solution of the system via Newton’s method. 

Given an initial estimate for the solution F % , a 0 , 6° we use Newton’s method to improve that 
estimate. The iterative process may be written as eq. (8), where AV, A / and J are defined in 
eqs. (9), (10), and (11). The Jacobian is given in terms of the dynamic pressure q , the incremental 
aerodynamic force vector, qR , the aerodynamic influence coefficient matrix, qA and the flexibil- 
ity matrix S. The incremental aerodynamic force vector is defined such that its component gr* 
represents the change in F a i due to a unit change in a, and the aerodynamic influence coefficient 
matrix, is defined such that its component ga,y represents the change in F a i due to unit change in 
0j. Similarly, the flexibility matrix, is such that is the change in 6 ,• due to a unit change in F a j. 


• solution by Newton’s method 


where 


JAY = A f 
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SOLUTION PROCEDURE (continued) 

Partial solution of equation (8) yields the following three equations for the increments AO, A a 
and AF a , shown below as eqs. (12), (13) and (14). We start with a rigid wing approximation and 
execute a single Newton iteration to approximate the flexible wing response. 


(I - qSA x )A0 = SBAh + 
Aa 

AF a 


SR 


A/a - N T Afi - qN T AA9 
q N T R 

Afi + qRAa + qAAO 


where 


B — I — 


rn t 

N T R 


A X = AB 


(12) 

(13) 

(14) 


(15) 

(16) 


• initial conditions 


F — 

1 ar — 


Oi r = 


/i(p,0,0) + qa r R 

jnjV — N^fi (p, 0, 0) 

qN T R 


(17) 

(18) 
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SENSITIVITY ANALYSIS (modular approach) 


As j ust stated, it is common practice to follow this procedure and use a single Newton’s 
iteration in the analysis of a flexible wing. Then for a design problem, where derivatives with respect 
to a design parameter p are required, equations (12), (13) and (14) are differentiated with respect 
to p (e.g., Ref. 2). This approach requires the calculation of derivatives of the matrices A and S 
which can be very costly. Here, instead, we follow Ref. 3 and differentiate equations (3), (4) and 
(7) with respect to p to obtain eq. (19), where a prime denotes differentiation with respect to p and 
where Y' and /' are defined in eqs. (20) and (21). The Jacobian J appearing in equation (19) is 
the identical matrix utilized in the analysis in equation (11). Equation (19) can be partially solved 
to yield the sensitivities shown in eqs. (22), (23) and (24). 

This approach does not require any derivatives of A and S but only partial derivatives of /i, 
fi and fz. For example, f[ denotes the derivative of F a with respect to a design variable when a 
and 0 are fixed. 


• modular sensitivity 

JY' = /' 

where 

Y‘ = |F' a' 6’f 

/' = [/; s'2 Kf 

• partial solution 

CD 

(I - qSA*)t>' = SBfi + + f' 3 

, fi - N T f[ - qNTM 
“ qN T R 

F' = / ; + qRa' + qAO' 


( 19 ) 

( 20 ) 
( 21 ) 

( 22 ) 

( 23 ) 

( 24 ) 
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SENSITIVITY ANALYSIS (traditional approach) 


By contrast, the more traditional approach (e.g., Ref. 2) to the derivative calculation is 
obtained by differentiating the aeroelastic analysis equations, such as eqs. (12) to (14) with respect 
to p as shown in eq. (25). This complicated expression can be shown to be equivalent to eq. (22). 
However, the traditional approach which employs eq. (25) requires the expensive calculation of the 
derivatives of the aerodynamic influence coefficient matrix, A' and the derivatives of the flexibility 
matrix S'. 


(/ - qSA x )A0' = qS'A x A6 + qSA'BAO + qSAB'AS 


+ S' B Aft + SB' Ah + SBA }[ 

+ WR Af2 + S ^WR ^ Af2 + Wr 


+ A/3 


( 25 ) 
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AEROELASTIC DIVERGENCE 


The aeroelastic divergence instability is calculated at a fixed angle of attack, because it is 
assumed that the pilot does not react fast enough to change the angle of attack as the wing 
diverges. The instability is characterized by a homogeneous solution to eq. (8), that is given in eq. 
(26). Equation (26) is an eigenvalue problem for q. The lowest eigenvalue is the divergence dynamic 
pressure qD • Equation (26) can be reduced to a standard linear eigenproblem by substituting for 
A0 in terms of A F a to obtain eq. (27). We denote the solution of eq. (26) as [F a D,0D] T and note 
that the same eigenvalue problem has also a left eigenvector [F a L^L] T • 


• eigenvalue problem 



(AS - ij)AF a = 0 (27) 

• right and left eigenvectors 


I -qoA 
-S I 

r kd] 
{ to J 

[ = ° 

(28) 

1 l 

—Qd^ 

I 

= 0 

(29) 
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DIVERGENCE SENSITIVITY 


To find the derivative of the divergence dynamic pressure qo with respect to a design parameter 
p, we differentiate eq. (28) at q = q D with respect to p and obtain eq. (30). We premultiply eq. 
(28) by the left eigenvector, defined by eq. (29) and obtain eqs. (31) and (32). Equation 

(32) contains derivatives of A and S with respect to p which we have managed to avoid before. 
However, the corresponding terms can be simplified. Using the definition of S , eq. (11), we note 
the relationship in eq. (33). 

To see how S’F a D can be calculated without obtaining S' consider a more generic case. Let / 
be a function of a vector X , and let D be another vector. Let Xo be a particular choice for X , then 
eq. (34) provides us with a way of calculating the product df /dX(Xo) times D without calculating 
the individual components of df/dX. Therefore, to calculate S'F a n we start by calculating the 
derivative of / 3 to a perturbation in F a in the form of F a o (because we use linear structural analysis 
this is the response of the structure to F a o)- Then we calculate the derivative of this response with 
respect to p assuming that F a D is fixed. The term A'$o in eq. (32) is treated in a similar way. 


to find q' D 
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• to efficiently find S' F a D 


S'F aD 


- fL (* 1 l\ F 

dp U Fa) FaD 


with 


(30) 


(31) 


(32) 


(33) 


|L(Xo)D = for i 7 l/(X 0 + tD) - /(*„)] = ton f/(X 0 + <D) (34) 

c/A e €— ae 
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APPROXIMATE OPTIMIZATION PROBLEM 


The optimization problem addressed in this paper is to minimize the structural weight W of 
the wing subject to aerodynamic, performance and structural constraints. It can be written as eq. 
(35), where gi, gi and <73 denote aerodynamic, performance, and structural constraints, respectively. 
The vector of circulation strengths T is calculated from eq. (1) and the nodal displacement vector, 
U, is calculated from eq. (5). 

Even with the more efficient sensitivity analysis, a fully coupled structural-aerodynamic anal- 
ysis and sensitivity is quite expensive. Thus, it is not feasible to optimize the design problem by 
directly connecting an optimization algorithm with the analysis procedure. Instead, a sequential 
approximate optimization algorithm is considered to be the best approach (e.g., Ref. 7). This 
approach replaces the original objective function and constraints with approximations based upon 
nominal values and derivatives at an initial point. Move limits are used to prevent the design from 
moving outside the bound of validity of the approximations. 

The approximate optimization problem is based on a linear approximation of the aerodynamic 
and structural constraints about a candidate design point po- That is, the approximate constraints 
0io and gza. are given in eq. (36), where A p = p— po» The performance constraints are typically quite 
nonlinear and inexpensive to calculate, so they are calculated exactly from the linear approximation 
to the aerodynamic solution. The approximate optimization problem is given then in eq. (37), 
where E represents a vector of move limits imposed to guarantee the quality of the approximation. 


• optimization problem 

minimize W(p) such that gi(T,p) > 0 

g 2 {r,o>,p) > 0 
gs{u, P )> 0 

• approximate constraints 

9ia{p) = 9i{Po) + g'i{Po)&P 

93a{p) = 93 {Po) + 03 fao)Ap 

• approximate optimization problem 

minimize W(p) such that g\ a (p) > 0 

S2(r„,a„,p) > 0 
Q'iaiv) > 0 
||Ap|| < E 


( 35 ) 


( 36 ) 


( 37 ) 
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APPROXIMATE OPTIMIZATION PROCEDURE 


The approximate optimization problem is solved sequentially as shown in the flowchart below, 
until the change in the design is smaller than a specified tolerance or the improvement in the 
objective functionis smaller than another tolerance. After anoptimum is found, a new approximation 
is constructed there, and the process is repeated until convergence is achieved. The optimizer used 
is the NEWSUMT-A program, Ref. 8, which is based on an extended interior penalty function 
procedure, and allows for various levels of constraint and objective function approximations. 
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SENSITIVITY TIMING COMPARISONS 


This figure presents a comparison of derivatives of vertical displacements and divergence dynamic 
pressure with respect to one structural design variable using the modular approach and the direct 
approach. We see that the values are very close. 

For structural design variables, the modular approach is also shown to save 32% in CPU time. 
Larger savings are anticipated for aerodynamic variables, because these entail the more costly 
calculation of the aerodynamic influence coefficient matrix and its derivatives. 


O Sensitivity comparison of vertical displacements to skin thickness 


modular approach 

direct approach 

9.2006574E-06 

9.2006600E-06 

-7.9199166E-06 

-7.9186167E-06 

-4.0682905E-06 

-4.0674891E-06 

-4.4964883E-06 

-4.4956877E-06 


O Sensitivity comparison of divergence dynamic pressure to skin thickness 


modular approach 

direct approach 

1.7295838E-03 

1.7340794E-03 


O CPU comparison for 1 design variable 

• modular approach : 6.58 sec. 

• direct approach : 9.59 sec. 




CONCLUDING REMARKS 


This paper focused on the processes of simultaneous aerodynamic and structural wing design 
as a prototype for design integration. The research concentrated on the major difficulty associated 
with multidisciplinary design optimization processes, their enormous computational costs. Methods 
were presented for reducing this computational burden through the development of efficient methods 
for cross-sensitivity calculations and the implementation of approximate optimization procedures. 
Utilizing a modular sensitivity analysis approach, we showed that the sensitivities can be computed 
without the expensive calculation of the derivatives of the aerodynamic influence coefficient matrix, 
and the derivatives of the structural flexibility matrix. The same process was used to efficiently 
evaluate the sensitivities of the wing divergence constraint, which should be particularly useful, not 
only in problems of complete integrated aircraft design, but also in aeroelastic tailoring applications. 


O Modular approach applied to integrated design 
O Improved divergence sensitivity 
O Computational efficiency of the modular approach 
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